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Inert Doublet Model with a 125 GeV Higgs 

M. Krawczyk, D. Sokolowska, B. Swiezewska 
University of Warsaw, Faculty of Physics 

A 125 GeV Higgs-like particle discovered at the LHC in 2012 has properties expected for it in the 
Standard Model (SM), with a possible enhancement in the two-photon channel. Such SM-like Higgs 
scenario can be realized within the Inert Doublet Model (IDM) - a version of the Two Higgs Doublet 
Model with an exact discrete D (^2-type) symmetry. In this model one SU(2) doublet plays the 
role of the SM Higgs doublet with one SM-like Higgs boson. The second doublet has no vacuum 
expectation value and does not interact with fermions. Among four scalars constituting this D-odd 
doublet the lightest one is stable, being if neutral a good DM candidate with the right relic density. 
In this paper an analysis of the two-photon Higgs decay rate in IDM, respecting theoretical and 
5^ other experimental constraints, is presented. The enhancement in h — > 77 is possible only if invisible 

Oh channels are closed, with the enhancement _R 77 >1.2 for masses of DM and charged scalars below 

■^T 154 GeV. The temperature evolution of the Universe to the Inert phase (described by the IDM) in 

T 2 approximation is presented and all possible sequences of phase transitions in one, two and three 
steps are analyzed. Going beyond this approximation by using an effective potential approach with 
one-loop T = Weinberg-Coleman term and temperature dependent effective potential at T 7^ we 
address the question, whether the strong first-order phase transition needed for baryogenesis can be 
realized in IDM. A region with a good DM relic density and a strong first-order phase transition is 
established, and discussed in the light of the XENON-100 data. 
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I. INTRODUCTION 



The Theory of Matter refers to the basic concepts like the quantum field theory framework, gauge symmetry 
►^J and the assignments of the lightest fermions. The idea of Spontaneous Symmetry Breaking (SSB) also has such 
r — status, although a particular realization of the SSB is not known and various models are considered. Standard 
iy~j Model (SM) is based on one SU(2) doublet of spin-zero fields. Standard extensions of the SM are models with 
£■ — two doublets (2HDM), which can be realized in MSSM. Similarly, models with more doublets, as well as with 
singlets, can be considered, all leading to p = 1 at the tree-level. In 2HDM various versions of the Yukawa 
■^1- interactions are possible. 

There are five scalars in 2HDM, among them neutral SM-like Higgs particles. The new particle, discovered 
in 2012 at the LHC, strongly resembles the Higgs boson from the SM. Here we consider a particular version 
of 2HDM - the Inert Doublet Model (IDM), which is very similar to the SM, in the sense that only one SU(2) 
doublet is involved in the SSB and there exists only one Higgs particle with the tree-level properties like the 
Higgs boson in the SM, H SM . The properties of the second doublet are quite different: it is not involved in the 
SSB and does not interact with fermions. It contains four dark scalars, which have limited interactions with 
the SM particles and the lightest of them is stable, thus, if neutral, being a good candidate for the Dark Matter 
(DM) particle. 

II. D-SYMMETRIC 2HDM 

The real content of the theory is determined by the symmetry properties of the Lagrangian as well as of 
the vacuum state. Here we assume that the potential is ^-symmetric with respect to the transformation 
4>s —* 4>Si 4>d —*>■ —4>d, which we call below the D symmetry. This symmetry leads to the CP conservation in 
the model. 

The D-symmetric potential has the following form: 



V = -- 
v 2 






with all parameters real. We take A5 < without loss of generality pQ . 

Various extrema (which can be local or global minima or the saddle points) can be realized in this potential. 
The vacuum expectation values (vevs) are as follows: 



1 r°v u^=±( u 



vuv """"V5W vs > VD > ueR - (1) 
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Neutral vacua are realized for u = 0. There are four types of neutral extrema (so also vacua) that have different 
symmetry properties: 

(i) The EWs case with u — vr> — vg = corresponds to the EW symmetry. 

(ii) Mixed vacuum (M) with u = 0,ug ^ 0,vd ^ 0. There exist two charged Higgs particles H ± , a pseu- 
doscalar Higgs A and two CP-even Higgses h and H, both of them could be SM-likc. 

(iii) Inert vacuum (Ii) with u = vrj = Q,vg ^ is the only state that conserves the D-parity and assures the 
existence of a stable scalar particle - the DM candidate. 

(iv) Inertlike vacuum (I2) with u = v$ = 0, Vd 7^ 0, which spontaneously violates the D symmetry by Vd =/= 0. 

If u 7^ then the charged vacuum (CB) is realized. This leads to the XJ(1)qed symmetry breaking and the 
appearance of a massive photon. Such a case is not realized in the nature today. 

Mixed and charged minima cannot exist for the same values of the parameters of the potential V. On the 
other hand, minima of the inert-type (I\ or I2) can overlap one another and M or CB in the parameter space, 
e.g. in the (A^As) plane [2]. 

A stable vacuum exists only if 



Ai>0, A 2 >0, A 3 + \/Xi\ 2 > 0, A345 + v%M >0 (A345 = A3 + A4 + A5), 



(2) 



so that 



R = 



A345 



/A1VA2 



> -1. 



Extremum with the lowest energy corresponds to the global minimum, i.e. the vacuum. The vevs (TTj) are related 
to the parameters of the potential through the extremum conditions and so the various values of vs,vd,u can 
be represented on the phase diagram (/ii,^), where 



Mi 



ran 

%/A7 



M2 



/As 



Different regions of this parameter space correspond to the different vacua. There are three regimes of parameter 
R which correspond to very different phase patterns shown in Fig. [T] There is a possibility of unique coexistence 
of two inert-type minima, as shown in Fig. |la| 
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FIG. 1: Phase diagrams for the neutral vacua in the D— symmetric potential: regions of /ii, /i 2 where various 

neutral minima EWs (white), 1^ (yellow/light gray), I 2 (blue/dark gray), M (green/medium gray) can be 

realized are shown. The hatched region in the panel (a) corresponds to the coexistence of I\ , Ii minima. 



Analogous phase diagrams can be considered for EWs, Ii,l2, CB with R — >• R3 = X^/^/XiXi and M — > CB. 
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III. IDM 

The Inert Doublet Model (IDM) is realized if I\ is the vacuum state. Here only one doublet (4>s) is involved 
in the SSB and there is only one SM-like Higgs boson h - we assume that its mass is equal to 125 GeV. The 
doublet (f>rj has vev = and Yukawa interaction is as in Model I, since (j>rj does not interact with fermions (it is 
inert). The D symmetry is exact here and only the second doublet is D-odd. It contains four scalars H, A, H^, 
called also dark scalars (D-scalars) . The lightest neutral scalar H (or A) may play a role of the DM. We take 
DM = H (so A4 + A 5 < 0). 

Masses of the scalar particles read: 

Ml = W = m\ x , Ml ± = hv " ~ m ™ , Mi = A4 ± + ^V, Ml=Ml ± + ^±^v\ 

with v = 246 GeV. 

These masses and the physical couplings can be used to describe the properties of the IDM. A345 is related 
to triple and quartic couplings between the SM-like Higgs h and the DM candidate H . A 2 gives the quartic 
DM self-couplings, while A3 describes the Higgs particle interaction with charged scalars. These parameters 
are subject to various theoretical and experimental constraints (see e.g. |2Tll7|): 

a. Perturbative unitarity constraints on self- couplings Parameters A of V are constrained by the condition: 
I A,; I < 8tt |13ffT5] . where Aj are the eigenvalues of the high-energy scattering matrix of the scalar sector. This 
leads to the following limits: 

A^ x = 8.38, A 3 e (-6.05, 16.53), A 345 e (-8.10, 12.38). 

b. Constraints from existence of I\ vacuum Existence of the normal (Mixed) vacuum is equivalent to having 
positive scalars' masses squared (condition for the minimum). It is not enough in the case of the Inert vacuum, 
as a coexistence of Inert and Inertlike minima is possible (see Fig. [Taj . For the Inert state to be the global 
minimum in addition the following condition has to be fulfilled |1 : 

m\ x m\ 2 



/AT 



> ^- (3) 



Knowing the Higgs boson mass Mh — 125 GeV and using the unitarity limit A™ ax = 8.38 one can derive the 
following limit on m\ 2 , that was not considered before |15j : 

m§ 2 < 9- 10 4 GeV 2 . 

c. Electroweak precision tests EWPT strongly constrain physics beyond SM. Values of S and T parameters 
are demanded to lie within 2a ellipses in the (S, T) plane, with the following central values jTH]: S = 0.03 ±0.09, 
T = 0.07 ± 0.08, with correlation equal to 87%. 

d. LEP II limits The LEP II analysis excludes the region of masses where simultaneously: Mh < 80 GeV, 
M A < 100 GeV and M A - M H > 8 GeV. For M A - M H < 8 GeV the LEP I limit M H + M A > M z applies 

01 El. 

e. A2 parameter Neither there exist, nor are expected in the near future, experimental limits on quartic 
selfcoupling for dark scalars A2. However, as we pointed out in [19, 20J there are some indirect constraints for 
A2 that come from its relation to A345 through the vacuum stability constraints (|2| and existence of I\ vacuum 

©• 

IV. RELIC DENSITY CONSTRAINTS 

IDM provides a good DM candidate (H) in agreement with the data on relic density QrjMh 2 , ^dm^ 2 = 
0.1126 ± 0.0036 P1J, in three regions of M H [5HLH HH HOI 121 122 : 

• light DM particles with mass below 10 GeV. Other dark particles are much heavier than DM and so this 
scenario mimics the behavior of the singlet DM model. 

• medium mass regime of 50 — 150 GeV with two distinctive regions: with and without coannihilation of H 
with the neutral D-odd particle A. 
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• heavy DM of mass larger than 500 GeV. In this region all dark particles have almost degenerate masses 
and coannihilation processes between all dark particles are crucial for the agreement with the measured 
DM relic density. 

The relic density data can be used to constrain A345 coupling for chosen values of masses of H and other scalars 
IB]. In Fig. 0the QoMh 2 as a function of A345 for various values of the DM mass is shown. Horizontal lines 
correspond to the 3a limit 0.1018 < Q DM h 2 < 0.1234. 

For both cases, i.e. (a) with and (b) without coannihilation of H and A, shown respectively in Fig. 2a and 



Fig. 2b the allowed region of A345 is symmetric around zero for masses Mh ;$ 72 GeV. Usually small values of 
A345 are excluded due to a non efficient DM annihilation. In the case (a) allowed values of A345 are smaller than 
in the case (b), as the process HH — > // with the cross-section a ~ Ag 45 does not have to be that efficient to 
provide the proper relic density value. As mass increases, the region of proper relic density shifts towards the 
negative values of A345 , which is due to opening of the annihilation channels into the gauge bosons final state 
and interference of processes HH — > h — > VV and HH — > VV. 




-0.2 -0.1 0.0 0.1 

(a) M A = M H + 8 GeV 




-0.1 0.0 0.1 

(b) M A = M H + 50 GeV 



FIG. 2: (A345, UrjMh 2 ) for various values of Mjj with (a) and without (b) coannihilation. Horizontal lines 

denote 3<r WMAP limit: 0.1018 < n DM h 2 < 0.1234. For tt DM h 2 < 0.1018 H from IDM does not constitute 

the whole DM in the Universe. Region of Vtr>Mh 2 > 0.1234 is excluded. Numbers on the curves denote the 

DM mass M H (GeV). We set M h = 125 GeV and M H ± = M H + 50 GeV. 

As it will be shown in the next section, in the IDM the enhancement of the 77 decay rate for the 125 GeV 
Higgs is possible for A3 < 0. If H is the lightest dark particle (DM candidate) then also A4 + A5 < and so 
we get A345 < for such case. For a recent analysis of the DM properties in the IDM and its relation to the 
two-photon decay rate of the Higgs boson see also Ref. |24) . 

It is worth noticing that the relic density as all other experimental data does not limit the A2 parameter. 



THE TWO-PHOTON DECAY RATE OF THE HIGGS BOSON 



The two-photon decay rate of the Higgs boson in the IDM reads: 

a(pp -*• ft -S- 77 ) IDM Br(/i 



R-, 



77) 



IDM 



o~(pp — > h — > 77) 



SM 



Br(/i -> 77) 



SM 



where the fact that the main production channel is the gluon fusion and that the Higgs particle in IDM is 
SM-like (in the sense defined before) was taken into account. In the IDM there are two possible sources of 
deviation from i? 77 = 1: 

• A contribution from the charged scalar loop to the partial decay width T(h —t 77) IDM 

MDM _ G F a 2 Ml 



T(h -► 77) 1 



\M 



SM 



SM 



IDM I 



128V27T 3 
where M SM is the amplitude present in the SM and 5M IDM = 2M f,% +m ' 2 v4 



LA./: 



Mr. 



with 2M I 



H±~ 



A3W 2 (the definition of the function Aq can be found in Ref. 



]). The charged scalar contribution can 
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interfere either constructively or destructively with the SM contribution, this way increasing or decreasing 
the decay rate. 

• The invisible decays: h — \ HH and h —} AA can augment the total decay width T lDM (h) with respect to 
the SM case. 

Performing a scan of the parameter space [30[ , taking into account the theoretical and experimental constraints 



presented in Sec. Ill we found the regions where i? 77 >l, with the maximal value of i? 77 around 3.4. Fig. 3a 
shows values of R 11 as a function of Mjj - it is clear that i? 77 > 1 is not possible for Mh < Mh/2 ~ 62.5 GevT 
It means that if the invisible channels are open, they suppress the charged scalar loop effect and enhancing i? 77 
is impossible. 
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(b) Region allowed by the theoretical and 

experimental constraints in the (m| 2 : ^ff±) plane. 

The curves correspond to the fixed values of i? 77 

(for the invisible channels closed). 



FIG. 3: Results on i? 77 for IDM, points with i? 77 < 1 (i? 77 > 1) are displayed in dark green/gray (light 

green/gray). 



Fig. 3b shows the viable parameter space in the (rn\ 2 , M H ±) plane together with the curves corresponding to 



constant values of i? 77 (calculated for the case with invisible decay channels closed) . It can be seen that than 
the i? 77 enhancement is possible only for: 

m 2 2 <-9.8-10 3 GeV 2 , 

which is equivalent to the bound A3 < 0, and consequently A345 < 0, see also Ref. [29] . We found that i? 77 > 1 
can be achieved for any value of M H ± . However, if larger value of i? 77 is demanded, then allowed values of M H ± 
are constrained. For example, for i? 77 > 1.2 we get the following bounds on M H ± and Mh (as Mh < M H ±)- 

62.5 GeV < M H < 154 GeV, 
70 GeV < M H ± < 154 GeV. 

VI. EVOLUTION OF THE UNIVERSE TOWARDS THE INERT PHASE 

We assume that today we live in the I\ vacuum, however in the past the vacuum could have been different. 
We first consider the simplest description of the evolution of the Universe with the thermal T 2 corrections. In 
this approximation, valid for high temperatures T 3> M , A parameters are fixed and only mass parameters vary 
with temperature as follows: 



AT) 



c t T 2 (* = 1,2), 
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where 



c\ = 



3A X + 2A 3 + A 4 3g 2 



9' 2 , 9 2 t +9t 



(» 



8 



C2 = 



3A 2 + 2A 3 + A 4 3.g 2 + g' 



12 



(i 



Here g, c/' are the EW gauge couplings, while gt,<?b are the SM Yukawa couplings |TJ |3T]. From the positivity 
conditions and the positivity of the vevs squared it follows that c\ + c 2 > 0. 

The quadratic coefficients of V change with temperature T and so the ground state of the potential may 
change. There exist three types of evolution of the Universe from the EW symmetric state to the I\ phase in 
one, two or three steps: EWs — > I\, EWs — > 1% — > I\ or EWs — > 1% — > M — > I\. In general, in the considered 
T 2 approximation, phase transitions are of the second order with the exception of the 1% — > l\ transition in the 
two-stage evolution. This is a transition between two degenerate minima I 2 and ii and it can be realized only 
if R > 1. Notice, that the DM exists only in the D-symmetric state I\ and so it may appear later during the 
evolution, if the Universe goes through the two- or three-stage sequences. 

If R < only one sequence that corresponds to the restoration of EW symmetry in the past is possible, 
namely EWs — > I\ (Fig. Ba|. It can be realized when R 11 > 1, which is suggested by the recent LHC data, as 



discussed above 
(Fig 



4b I 



In the other scenarios for R < in the initial state of the Universe there is no EW symmetry 
There is a possibility of a temporary restoration of the EW symmetric state (scenario Y) , in other 
cases the EW symmetric state never existed. 

For a certain parameter range there is a theoretical possibility of having a charged vacuum in the past 
(Fig. [4c| [T]. This type of evolution is equivalent to having IDM with the charged DM particle PQ. Current 
model-independent bounds require the charged DM to be heavier than 100 TeV [32_. Heavy charged scalar may 
appear in the model without breaking the perturbative unitarity conditions for very negative TO§ 2 , i.e. large 



|m| 2 |. However, the sequence Z+ requires that c 2 < and |c 2 |/ci > \m\ 2 \lm\ 1 ^ (10 5 -r 10 6 ). This contradicts 
the requirement c\ > — c 2 based on the positivity condition and invalidates the possibility of passing through 
the charge breaking phase during the evolution pQ. 



W(T) 




M T ) 



, I2 (T) 



l'i{T) 




W(T) 



(a) R < 0: the only sequence that leads (b) R < 0: possible sequences with 

to restoration of EW symmetry in the non-restoration of EW symmetry, 

past. 



(c) R3 < 0: transition through the 
charged vacuum. 



FIG. 4: Possible sequences of phase transitions for R < 0. 



VII. TEMPERATURE EVOLUTION BEYOND T 2 APPROXIMATION 



The T 2 approximation, used in the previous section, simplifies the discussion so that all phase sequences, 
that are in principle possible, can be recognized and analyzed in details. However, this approach has obvious 
limitations, since some sequences are realized at relatively low temperatures, where the T 2 approximation is 
not valid. The analysis of the evolution of the Universe to the Inert phase going beyond this approximation was 
performed in [!JT], where the effective potential approach was used under the assumption that only neutral vacua 
(phases) were realized in the past. It was found, as expected, that some of the second-order phase transitions 
change their character to the first-order. The analysis which focused on the strength of the first-order phase 
transition needed for baryogenesis was performed, based on the one-loop effective potential at T = (Coleman- 
Weinberg term) and temperature dependent effective potential at T ^ (including sum of the ring diagrams) 



Toyama International Workshop on Higgs as a Probe of New Physics 2013, 13-16, February, 2013 

[33]. Of course, IDM as a model with the CP-conserving potential can not serve alone as a realistic model 
for baryogenesis and an additional source of CP violation has to be included. Nevertheless, one can examine 
whether the model parameters, for which the right density of relic particles is predicted, are compatible with 
the strong first-order phase transition that could allow the electroweak baryogenesis. 

The strong first-order transition is realized if for the critical temperature Tew we have v(Tew)/Tew > 1- 
We found that this condition is fulfilled in the M/j=125 GeV case with medium DM mass Mjj ~ 60 GeV and 
relatively heavy dark scalars A and H ± , with masses ~ 300 — 400 GeV (we took them to be degenerate for 
simplicity). The important parameter is A345 - the coupling of the Higgs particle h to the DM (hHH), while the 
parameter A2 does not play any role in the analysis (we fix it to be equal to 0.2). For the mass of the DM equal 
to 65 GeV we obtained allowed mass region of heavier dark scalars between 275 and 380 GeV (Ma = Mh± ) , 
and in order to be in agreement with XENON- 100 data | A345 1 should be smaller than 0.1. This is in agreement 
with similar analyses |34l [33], see also papers for 2HDM (not IDM) [36 39]. It is worth noticing that these 
regions of parameters correspond to the one-stage evolution of the Universe, and the negative A345 corresponds 
to the enhanced h — ¥ 77 decay rate. 

VIII. SUMMARY AND OUTLOOK 

2HDM is a great laboratory of BSM. A distinguished version of 2HDM is the Inert Doublet Model with 
the SM-like Higgs sector and the dark sector with a good DM candidate. Recent LHC data with the SM-like 
Higgs signal, with a possible enhancement in the 77 channel, put strong constraints on masses of DM and other 
dark scalars, as well as the self-couplings, especially A345. A striking conclusion is that an enhanced h — > 77 
rate excludes a light DM particle and sheds light on the evolution of the Universe, suggesting only one phase 
transition from the EWs to the today's Inert phase with the right DM relict density. 

The precise measurements of the SM-like Higgs properties are crucial not only to establish details of the EWSB 
mechanism, but also to understand the composition of the Dark Matter. The future linear collider in the e + e~ 
and 77, ej options may indeed lead to further progress in these directions |4"0"l |4"T] . It is known that for example 
the partial width T(Hsm — > 77) can be measured at linear collider with ~ 2% precision, by combining the 2% 
measurement of the Higgs particle production in 77 collision decaying into 66, T(Hsm — > "fj)Br(HsM — ► 66), 
with the measurement of the Br(HsM ~^ 66) to be performed with 1% precision at the e + e~ option of ILC at 
500 GeV [42 . Some analysis on IDM at the ILC has just appeared [43 . 
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